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Abstract 

We have constructed a sequence of solutions of the Helmholtz equation forming an 
orthogonal sequence on a given surface. Coefficients of these functions depend on an 
explicit algebraic formulae from the coefficient of the surface. Moreover, for exterior 
Helmholtz equation we have constructed an explicit normal derivative of the Dirichlet 
Green function. In the same way the Dirichlet-to-Neumann operator is constructed. 
We proved that normalized coefficients are uniformly bounded from zero. 

Keywords: explicit solution, Helmholtz exterior problem, Green function, Dirichlet-to- 
Neumann operator. 
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1 Introduction 

Consider Q C M 3 with Lipschitz boundary dQ and k > 0. The scattered field is given by 
Helmholtz equation and radiation condition 

R 3 \n, (i.i) 

R^oo, (1.2) 
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with Dirichlet boundary conditions, 

*(r) = uo(r), r g u e (7(90). (1.3) 

For example, in [1] is proved the existence and uniqueness of the solution of fll.ll) - fll.3p . A 
function \&(r) which satisfy mentioned conditions has asymptotics 

z/cjT*! / "j^ *\ 

%) = tt/(?)+° n . r -> oo, g = r/|r|GS' 2 , (1.4) 
\r\ \l r l/ 

where function f(9, <p) = f(9, <p, k, Uq) is called scattering amplitude and observable 

°t= I \f(q)\ 2 da(q) 
Js 2 

is called Total Cross Section, a is a square element of the unit sphere. 
A very important particular case of the boundary condition is 

UQ = e ifc <^°>, r G dtt, 

which is the scattering of a plane wave with incident angle 9$ G S 2 . The total momentum 
transmitted to the obstacle is given by observable called Transport Cross Section (in a large 
volume normalization) 

R = [ (l-<q,6 >)\f(q)\ 2 da(q). 
Js 2 

Unfortunately, analytical expressions of these observables for certain k > exist only 
for few bodies of simple shapes (see [TBI HE])- Moreover, the scattering happens not only 
by plane or spherical wave, but it could be caused, for example, by arbitrary secondary 
radiation. 

From another point of view there exists some numerical methods for direct scattering 
calculation. One of them is based on a numerical solution of integral equation (see |15j). 
Another method, developed by A. Ramm and S. Gutman in (Sj-fO], allow to construct the 
Green function and therefore to obtain solutions for arbitrary boundary condition uq G 
C(dQ). The ground analytical achievement by A. Ramm and S. Gutman is the so called 
Modified Rayleigh Conjecture. In particularly, it follows that functions Yi m (9, (p)hi(k\r\)\gn 
(spherical harmonics and spherical Hankel functions (see fll.7jl ) correspondingly) form a basis 
in the space L 2 (dfl, dS). 

The aim of current work is to express explicitly coefficients of physical observables and 
also for the normal derivative of the Dirichlet-Green function. We found a constant uniform 
lower bound for normalization coefficients (denominators) and we prove convergence of all 
produced series. 

Greens function of the Laplacian in Q' is given by 

A r G(r,t) + k 2 G(r,t) = -5(r -t), r,teQ', (7 = 0, rGffi, 
lim / 
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As it is known (evidence from Green formula) 

j ^p^- Uo (t)dS(t), reTV 
Jan on* 

In this paper we shall express kernel 9G d ^f explicitly through coefficients of the surface 
dfl and also as measures C, o~t, R and {Ai m }, where A m ,0 < \m\ < I are coefficients of 
representation 

oo 

;=o \ m \<l 

where Y m i(9,ip) are spherical harmonics. 

Let N = N U {0}. Let £ be the set of indexes 

£ = {(l,m) :I,meN, \m\ < 1} 

and let Z G £ with / = (1(1), 1(2)). We also we use the notation I = (1(1), —1(2)) and set the 

order: (I, m) > (p, r) < > I > p A [(/ = p) V (\m\ > \r\)] A [(/ = p) V (m > 0) V (m = —r)], so 

£ = {(0,0), (1, 0), (1, -1), (1, 1), (2, 0), (2, -1), (2, 1), (2, -2), (2, 2), (3, 0), ...}. Let o = (0, 0) 
and let operations + and — have the natural definition in £ x Z — > £ correspondingly to 
the introduced order. 



1.1 Surfaces with inverse radius- vector represented as finite com- 
bination of harmonics 

Let T be a subset of functions (multiindex) £ N which have finite support and let capacity 
be defined by 

\d\=J2d(l), Suppd = mzx{l:d(l)^0}. 
i 

Let ei G T be defined as ej(ra) = 5i m (evaluated 1 only when I = m). Also set 

C d _ W n d-Y[ n d{l) 



Jd= r r Y\(Y t (e, v )f l) ded v . 

Jo Jo , 



Theorem 1.1. Let a star shaped surface dQ be given as a set {r = r(9,ip) G 1R 3 ,# G 
[0,7r],</?G [0,27r)} where 

Wv?)l = v IfT^T (L6) 

with N > and where {ai m } are coefficients. Then 
1. Functions 

^™( r ) = ^2 c -nkY k (9,ip)h k{ i)(k\r\), n G £ 

k<n 
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satisfy ( U.ij ).( TOj) and their restrictions {^f n \dn, n £ £} /orm an orthonormal basis in 

L 2 {dn,d6d(p). 

Here 

c nn = 1/A„, n G N 

2 
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where coefficients h n j are defined from the well known representation for Hankel spherical 
functions J7P) /. 
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5. Moreover, exists numbers C{ = Ci(k,Q),i = 1,2 stzc/i iaai 



o < k < n, k,n G £. 
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Note 1. Define {A^, m, L E £} where A^ = for L < m and 

■A m ~ ^ ^ Cnm ^ j C n pUpj L ^ Til . 
L>n>m V^L n 

or, recursively 

A m = A m + —CL m ^ ^ CLpUp. 

p<L 

From Theorem 11.11 and using the orthogonality of basis {^nla^} we have 

C 



\Am — A m \ < 



m(l)! 
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-E 
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o(l), L — > oo 



Then 
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lay — 0~k\ < C 



p<L 



o(l), -L — > oo . 



For the case m = e , i.e. 6*0 = (0, 0, 1), we have 



E 



(n(l) + n(2))(n(l)-n(2)) p 
(2n(l) + l)(2n(l) - 1) ^^"^a. 
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\R — Ri\ < \\u — u p \\ = o(l), L 

p<L 

Formula (11.131) is a simple corollary of equality xP n (x) = ^^P n+ i(x) + -^^P n _i{x) f° r 
Legendre polynomials. 

Note 2. The Dirichlet-to- Neumann operator M (which correspond the boundary con- 
dition Mo to the normal derivative of the field u(r)) could be constructed in the following 

way: 

(Mu ) (r) = £) ( I' r ^ q M*«WM)M*p\ ^r(r). 



jo 
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1.2 Comparison with other explicit representations 

In spite of Theorem 11.11 be stated for star shaped bodies clearly all representations of the 
Green function will be the same for arbitrary, smooth enough, body. We just should put 

pit i*2t: 

u P = / u (6,ip)Y p (6,ip)h p{1) (k\r(d,ip)\)ddip, p e C, 
Jo Jo 

9ij= / 2 V i (^,^)F J -(0,^)/ ii(1) (A;|r(^,^)|)^ (1) (A;|r(^,^)|)^, i,j G C 
Jo Jo 

It is possible, like we mentioned above, functions {Y p (9, <f)h p ^(k\r(9, <f)\)} do form a basis 
for arbitrary, smooth enough body. 

Our representation holds for all values of k > 0, but exists other explicit representations 
for the Green's function which holds in a neighborhoods of points k = and k = oo. We 
consider of interest to cite the following 

Theorem 1.2. R. E. Kleinman /2l Th. 4- J] There exists a > such that when \k\ < a, 
the Green's function G(x,r) exists uniquely in Q' and is given explicitly by 

G(r,t) = -i-r—> +e^J2 Kn ^ U- 14 ) 



n=0 



where 



K(U ) = -2ik [ dvih f ^ f . [\h\Uo], 
Jn' l J i| a \H\ 



r e ~ik\to\+ik\to-t\ Q 

U = U Q (ti,t)= / dS(t )—— __Gf (ti,to). 

Here v is a volume element, function Go{r,ti) is the static Dirichlet- Green function. 

As we can see, all summands here are determined recursively as in the Theorem 11.11 
Also, is important to note, that the first approximation is a static Green function, which is 
supposed to be known. So, in the some sense the constructed solution is a pertubated static 
solution, that explain the quick, exponential, convergence of the approximation (I1.14p . 

In the high-frequency exists the following representation for the Green's function which 
holds for obstacles which satisfy non-trapping condition and single impact condition (see [12] 
for details) 

ik\r— 1\ ^ / * \ 

GM) = ~i^~t\ +elkS(r,t) £ (£) **.(r,t) + J*(r,f), (1.15) 

' ' m=0 ^ ' 

where function S(r, t) satisfy eikonal equation 

|V r S| 2 = l, 

and function z m (r,t) are defined recursively through differential (transfer) equations: 
2V r S ■ V r z m + (V r S)z m = —V r z m _i, m > 0, 2_i = 0. 
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If r and t don't lie on the tangential rays, we have the estimate 

\Ri(r,t)\ = oik 1 ' 1 ), k^oo. 

In that case, the good, exponential, convergence exists since solution at high values of k are 
close to the "classical" limit solution. 

ot representation fll.Sp holds for all values of k > 0, but we don't have any estimates 
about the rate of convergence. For the sphere, our representation is the well known solution 
(for ex. [19]) with very good convergence (as a 1/nl). Of course, good convergence properties 
should be kept in a some neighborhood of the sphere. Moreover, we are sure that it will be 
possible to do it with developing of local limit theorems for sign distributions, really if 1/r 
is a characteristic function of a random variable with a sign distribution then elements of a 
Gram matrix g^ are in the sense of the probability of sum of (i+j + 2) of such variables (see 
[TT] for integral theorem). By now it is out of our possibilities to describe the neighborhood 
of at least exponential convergence. From another point of view we want to remind that 
representation (11.81) . in contrast to representations (11.151 II. 14}) . converges for all values of 
k > 0. 

To finish our discussion we want to note that representation like (jl.8p can not be done 
for the 2dimensional model, since Bessel functions of integer index do not allow a finite 
polynomial representation. 



1.3 Polyhedrons and bodies which are revolutions of polylines 

Calculation of the Gram matrix G = {gij} for polyhedron is very natural and elements g^ 
need only N (i + j) 2 sums, where N is the number of sides. 
Let's set iV + 1 points on the segment [0, n] 

= 9 < 9 1 < . . . < 9 N ^ <9 N = n, 

and consider the set of bodies 1Z which contains arbitrary bodies of revolution of the polyline 
with vertexes in points {6^, i — 0, . . . , iV}. 

Theorem 1.3. Consider polyline r = ^\ ri(9)Ii(9), where Ii,i = 1, . . . , N is a characteristic 
function of the segment 6j\ and coefficients {a^, bi, Ci,i = 1, . . . , iV} determine parts of 
the polyline ri: 

airiCos9 + birisin9 = fi, i—l,...,N. (1-16) 
Then Gram matrix coefficients gij could be calculated as 

~ ( f \i+j+2 Ly i+j+2 u p u p 1 ij,l{i+3Y 

p=l W 1=0 

where 

r e p _ 

Tf lm = / Yi{B)Yj{9) cos 1 6 sm m+2 - 1 9 sin 9d9, I < m + 2; 

N i+j . / m _ ^ \ /m+2 

9ij = ^2 J^+2 I hinhj( m -n) , „ y m+2 1 I C m+2 a p b p + Tf j>lm 

p=l m =0 \n=0 ^ Jp> J \l=0 
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1.4 Explicit dependence of Gram matrix on frequency 

We can represent elements of the Gram matrix {gij,i,j G £} for arbitrary smooth enough 
surface as a polynomial of the inverse frequency k^ 1 , 

m=0 

where 

P? j = J2 h W h m( m -i)J Q dej^ dtp | r(g y) ]l; 2 , 0<m<i(l)+j(l). (1.18) 

2 Proofs 

Denote 

*„(r) = y w (0)^ n(1) (A;|r|), ne£, r£fl', 

where functions ^ n {f) satisfy ( 11. ip and ( 11.21) . Moreover, they have asymptotics at the infinity 
(see(0)) 

# w (r) ~ iy n (0)^-l |r|-oc. 
K |r| 

Let A be a linear operator acting from L 2 (dQ,dS) to L 2 (S 2 ,da) which corresponds u\qq to 
scattering amplitude f(q) (see (11.40 ). Here dS and do - are standard metrics on dQ and unit 
sphere S 2 . In jH 2.12,2.16], [31 2.10,2.13] it is proved that A is a bounded operator and in 
particulary it is proved that functions ^ n \an form a basis in L 2 (dQ,dS). So we have the 
transformation 

Let us prove that exists a constant C\ = C±(k, Q) such that 

K\<-%t, neC. (2.1) 
n{l)\ 

Consider sphere Sr := {\r\ = R} such that the body fl is strictly embedded in that sphere. 
Define also the A\ : L 2 (dQ, dS) — > L 2 (Sr, dap), where g?(Xr = i? 2 sin 9d9dip, as 

/" dG(x,y) 

(Aiu )(x) = / — ^ — u {y)dS{y), x e Sr. 
J on an 

Since function ^ has singularities only for x = y (for example, see [8]) operator Ax is 
bounded. Note that functions h n {k\r\) are constant on the Sr, 

4:7iR 2 \h n (kR)\ ■ \c n \ = |(AiWo,^n.)i 2 (s«)| < 

ll^iMo|U a (s«)ll y »IUa(SH) = 4tt.R 2 || A^qW^Sr) < ^irR 2 (I Ai || |K \\l 2 (8Q) ■ 
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So 



IcJ < 



l^i|IH M o| 



L 2 (sn) 



K(kR)\ 



taking into account the asymptotic (see [T71 pp. 358-364]) of h n {R) ~ jm^yc-i n ~ * °°> an d 
since R could be chosen to satisfy kR > 2, we obtain (12.11) (that corresponds to ( 11.101 in 
Theorem 1)). 

Let g(9, ip) be a density 



dS= [g(9, V )] 2 da = \r\ 



\ 



r 2 + 



d\r\ 



sin 2 9 + 



d\r\ 

dip 



d6d<p 



(2.2) 



and for simplicity we shall use the same notation da = d8dip for measure on the dQ. It is 
evident that if function / G L 2 (dQ, da) then / '/ ' g G L 2 (dQ, dS). 

Denote = ^nldn/g and construct an orthonormal basis ^ in L 2 (dQ,dS). We shall 



construct it in form: 



fc=o 

n-1 



(2.3) 



where 



Using 



we have 



Denote 



fc=o 



Tl-1 



i<ir-£i«,* 



di 2 



k=o 



1/2 



p=o 



ra— 1 fc 
fc=m p=o 



/An, m < n . 



* w (r) = 0§° = ^ c„ fc * fcl n G £, refl'. 

fc=o 

and let «o be a function in L 2 (dQ, da) and also u$/ g = Y^q Mn^J,- Then function 

oo 

u(r) = 5^u n $ w (r), rGR 3 \Q 

ra=o 

satisfy ffTll-fOl). 

Using (I2.3P we obtain the corresponding scattering frequencies (see ( II. 5p ). 



1 x ( n 

T 2j C n™ I 



vp=o 



where u p = {u /g, *°) L2 (dn,dS) = («o, *plen)^(Oi.«fc)- 

By construction J2p= Cnp u p — u n and that follows (11.111) . 
Let us prove now (11. . Set 



n-l 



\ n = \\F n \\ L2{ 9u, dS ), where F n = *° - §° fc )§° . 

fc=o 

By construction (AF n ,Y n ) = |, therefore ||AF n || > j, so 

||AF n || < \\A\\\\F n \\ L2idQ>dS) . 
Due to the representation of r(9, cp) and (12.2ft exists 5 = max^ |<7(0, | 2 < oo, so we have 

Xn=\\F n \\ M oa,da) > ^\\F n \\L 2 {an,dS) > (^H A|| A;)" 1 . 
Now, let us prove (ll.Sp . We write (see (II. 9p ) 

l^-l 2 = ^ E 



k 2 



p<n \ P<n m<l<n p<l 



^ oo 

fc2 E E 

Next, we calculate (\1>9, ^ 



ra=o m<n 



vp<n / \ P<ra m<l<n p<l 



Up 



2tt 



^ = (^,«5)= / d0/ ^(^^(^^^(fclr^^D^D^Ir^^l) 



'0 ^0 
f*7T /*2tt 



1 ft i(l)P 



^Y»l*J mm *„) E^ E 



p=0 v 1 17 / \p=0 v 1 17 



' l)m L dd L WMrtHW E jb^(EW*a: ) i,,, M „„ +2 



m=0 \ 2=0 

»(1)+J'(1) m _ ^ -27T 



|r(0 >V O|» 



^ (2) IT ^ E W^-o jT * f ^SS^ = (2 ' 4) 

m=0 2=0 u u ' ^ ' r / 1 

*(l)+i(l) 1 / m _ \ 

MP' E ^ E^Awm E cv/-+«»i. 

m=0 \2=0 / d:|d|=m+2,Suppd<(V,V) 

Theorem 11.11 is proved. 

Proof for Polyhedrons — Coefficients gr^. From ( 11. 16ft we have 

1 dj cos 6* + 6j sin 6* 

- = 7 , t = l,-..,N. 
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First of all note that 

N 



Develop gf jt 



9ij = 1t,9iP with^. := ( P %(r p (e),8)*j(r P (8),8)s™8de. 

p=l Jxp-i 



gf,= / Y^Y^h^krie^hjiklrKe^sinede 



i+j 1 m 



1^ 1^+2 1^ h m h j(m-i) / ^ 



o \r(6,tp)\ 



m+2 



m=0 Z=0 
i+j ^ m m+2 

1M Yl h HX)l h i(™-l) ( f \ m +2 Yl Cl m+2 al pK +2 ~ lT ij 
m=0 1=0 ^ Jp> 1=0 



Imi 



where 



2* = / ^(0)^(0) cos' 0sin m+2 ~' sin l<m + 2. 



Proof of the representation (I1.17P follows explicitly from (12. 4p . 
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